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Motivation 



• Further understand behavior of hot-wire measurements in icing conditions 

• Newer applications in ice crystals and mixed phase 



SEA Multi- 

Element 

Probe 




Particle impingement uniformity -> Rigby et al. (Aviation 2014) 
Temperature profile along wire Thermal model (this paper) 
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Outline 

• Describe Thermal Model 

- Governing equations 

- Model parameters 

• Results (model vs. experiment) 

- Dry (air flow only, no water) 

- Wet (i.e. cloud on with LWC or IWC) 

• Conclusions 
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Thermal Model Based on Energy Balance 
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Governing Equation 


d^Ts 1 dk 

i -I 

dx^ k dTs 


dTs 

dx 


P ka r n 8(JP . . i^p 

f Nu[Ts - Ta,o] - — [TV" " TslA + 


dgkAc 


kA, 


kAl 


1 d^y^Ql 

kAc dx 


Conduction Convection Radiation Heat Heat 

Generation Sink 

(Water /Ice) 


2"^ Order, Non-Linear ODE 
solved using MATLAB routine bvp4c 


Boundary Value Problem 
need two boundary conditions 
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Model Parameters 




• Convection correlations 

- (1 ) Sparrow et. al - 2004 

- (2) Generated from CFD 


• Probe operation 

- Maintains const, avg. temp 
(e.g. resistance) -140 °C 
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Results 



• Experimental results from 2012 NRC RATFac Tests 


• Dry Conditions (Po=13.5, 6.5 psia; 11=85,100,135 m/s) 

- Temperature profiles 

- Total power (experiment vs. model) 

• Effect of heat-transfer coefficient 

• Wet results (Pq = 13.5 psia, U = 85 m/s) 

- LWC Sweep 0-3 g/m^ 

- IWC Sweep 0-1 0 g/m^ 
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Results: Temperature Profiles - Dry 
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Experimental Data - Dry 
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Model vs. 


16-1 


083" - Wire 
Tinf = -5°C 
Nu = Sparrow 2004 
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Model vs. Experiment - Dry 
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Heating & Evaporation Model 
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Evaporation Efficiency 



Needed a way to estimate fraction of 
water that does not evaporate 


Defined a parameter called evaporation 
efficiency, r| 


evaporation potential 
incoming mass flux 


0<i;<l 


Use analogy of heat & mass transfer 

- Mass flux related to heat flux 


Evaporation area: 

- For 083,021 entire circumference 

- For HP, forward face of element 

Does not include bounce / splash 



Temperature profiles from CFD 
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Effect of LWC on Temperature Profiles 
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Wet Power - Effect of Conduction Losses 



‘wet ~ ‘total ‘dry 




qgenm 6.34 13.61 

qconv (W) -6.05 -6.05 

qcond (W) -0.28 -0.14 

qrad (W) -0.01 -0.01 

qwet (W) 0 -7.41 


• T raditional wet power calculation: 

P^et = 13.61 W - 6.34 W = 7.27 W 

However, energy needed to evaporate LWC was 7.41 W (1 .9% diff) 
Under measurement is due to conduction loss differences 
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IWC 
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Wet Power vs. IWC 




IWCt (g/m3) 
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Conclusions (1 of 2) 


• Developed steady-state hot-wire thermal model 

- Includes resistive heating, convection, axial conduction, radiation, and 
water/ice evaporation 

• Examined: 

- Temperature & power variation along the wire 

- Steady-state power 

• Model compared to SEA multi-wire probe data 

• For dry conditions: 

- Matched experiment to within: 

• 5.5% for 021 

• 9.2% for 083 

• 14% for HP 

- Max. conduction loss ~ 4% of total power for conditions examined 
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Conclusions (2 of 2) 


• Wet conditions: 

- Introduced “evaporation potential” to estimate % water evaporated 

• Needs validation 

- LWC: 

• Affected temperature profile of 021 most significantly; 

• In all cases, high evaporation potential, effect minimal 

• Conduction losses can be different dry vs. wet 

- For 021 at 3 g/m^ ^ 1 .9% difference in measured vs. actual 

- IWC(HPonly) 

• Model suggests a non-linear behavior of wet power and IWCt 

- Below 4 g/m^, linear relationship 

- Above 4 g/m^, non linear due to incomplete evaporation everywhere along wire 

• Limited available experimental data to see if trend is correct 

• Bouncing or splashing-type loss present in experiment complicate interpretation 


• Further examination & development of model planned 
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